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^ ; Abstract 

-)— > 

^ ' Conditions on the elastic stiffnesses of anisotropic crystals are de- 

rived such that circularly polarized longitudinal inhomogeneous plane 
'^ I waves with an isotropic slowness bivector may propagate for any given 

C ■ direction of the normal to the sagittal plane. Once this direction is 

Q ■ chosen, then the wave speed, the direction of propagation, and the 

direction of attenuation are expressed in terms of the mass density, 
the elastic stiffnesses, and the angle between the normal to the sagit- 
K^ ' tal plane and the normals (also called "optic axes") to the planes of 

!^ , central circular section of a certain ellipsoid. In the special case where 

this angle is zero, and in this special case only, such waves cannot 
^ ■ propagate. 

^' 

2 ■ 1 Introduction 

In classical linearized elasticity theory, a special role is played by infinitesi- 
rS ' mal longitudinal homogeneous plane waves. For such waves the amplitude 

ci . vector is parallel to the propagation direction n (say) so that all the parti- 

cles oscillate along that direction n and the motion is one-dimensional. Such 
waves may propagate in every direction in an isotropic compressible elas- 
tic material, but this is no longer the case for elastic anisotropic materials 
such as crystals. Possible directions of propagation of longitudinal homoge- 
neous plane waves, called "specific directions" by Borgnis |T] may be as few 
as three in an elastic anisotropic crystal. By assuming certain restrictions 
on the elastic constants, Hadamard [2] created a special model anisotropic 
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elastic material in which longitudinal plane waves may propagate in every 
direction. 

Here, we consider inhomogeneous plane waves. These waves are atten- 
uated in a direction different from the propagation direction. They may 
be described in terms of bivectors [3] - complex vectors - the amplitude 
bivector, A (say), and the slowness bivector, S (say), which may be written 
[1] S = NC, where the "directional bivector" C is written C = mm + in 
(m ■ fi = 0, m > 1, |m| = |n| = 1) and A^ is called the "complex scalar 
slowness" . Once the directional bivector C is prescribed, the slowness S and 
amplitude A are determined from the equations of motion. To prescribe C is 
equivalent to prescribing an ellipse with major semi- axis mm and minor semi- 
axis fi; this so-called "directional ellipse" for inhomogeneous plane waves is 
the analogue to the direction of propagation n for homogeneous plane waves 
[1]. The inhomogeneous plane wave is said to be "longitudinal" if A and S 
(and therefore also C) are "parallel": A A S = 0. What this means is that 
[3] the ellipses of A and of S (and C) are all "parallel" , being similar - same 
aspect ratio - and being similarly situated - parallel major axes and paral- 
lel minor axes. In particular we consider "circularly polarized longitudinal 
inhomogeneous plane waves" (CPLIPW). For such waves both C and A are 
isotropic, that is C ■ C = 0, A ■ A = 0; the ellipses corresponding to C and 
A are coplanar circles, the normal to the plane being a = m A fi. 

Here we seek to determine restrictions on the elastic constants such that 
for any choice of a, the normal to the plane of C and A (which is the plane of 
the motion), a CPLIPW may propagate. We call the corresponding materials 
"special" . 

Starting with the constitutive equation for a general anisotropic elastic 
crystal (which involves twenty-one independent elastic constants), we obtain 
necessary and sufficient conditions on the constants in order that CPLIPWs 
may propagate in every plane. It turns out that nine linear relations among 
the elastic constants must be satisfied so that the special model material 
has at most twelve independent elastic constants. For such materials we 
determine the general structure of the (symmetric) acoustical tensor. The 
complex slowness A^ of the CPLIPW is determined for all choices of isotropic 
C. Because the wave is circularly polarized, one eigenvalue of the acousti- 
cal tensor is double |1]. The remaining simple eigenvalue corresponds to a 
"transverse" inhomogeneous plane wave that is transverse in the sense that 
its amplitude bivector B (say) is "orthogonal" to C: B ■ C = 0, which means 
that the orthogonal projection of the ellipse of B upon the plane of C is a 



circle. 

The equation giving the complex slowness A^ for the CPLIPW is of pre- 
cisely the same structure as the equation giving the complex slowness of the 
transverse wave. Both are of the form C ■ 0C = iV~^, C ■ C = 0, where 
is a real symmetric tensor. An ellipsoid E (say) may be associated with 
[5]. It is seen that the slowness bivectors are obtained by first determining 
the central ellipsoidal section S (say) of the ellipsoid E by the central plane 
with normal a. Then m and n are chosen to lie along the principal axes of 
S. This fixes C, and then A^, and therefore S (= A^C) is determined. For 
choices of a along the normals to the planes of central circular sections of the 
ellipsoid E, there is no propagating wave. 

Finally, we briefly consider the possibility of having CPLIPWs propagat- 
ing in crystals of various classes: triclinic, monoclinic, etc. It is seen that 
CPLIPWs may not propagate in trigonal, tetragonal, and cubic crystals, nor 
in isotropic materials. It is seen that they may propagate in triclinic, mon- 
oclinic, orthorhombic, and hexagonal crystals, provided the linear relations 
among the elastic constants evoked above (or their specialization to those 
classes of symmetry) are satisfied. As kindly pointed out by a referee, we 
must ensure that the crystals are purely elastic so that mechanical fields are 
not coupled to electrical fields. Hence in what follows, we restrict our atten- 
tion to the following crystal classes: 1, 2/m, mmm, 3m, 4/m, 4/mmm, 6/m, 
6/mmm, 432, m3, m3m. 

2 Basic equations 

The constitutive equations relating stress components aij with strain compo- 
nents Cij for a homogeneous anisotropic elastic crystal are given by Hooke's 
law: 

aij = dijkiCki- (2.1) 

Here dijki the elastic constants, or stiffnesses, have the symmetries, 

dijki = djiki = d^iij, (2.2) 

so that there are at most 21 independent stiffnesses. Also, 

2eij := dui/dxj + duj/dxi, (2.3) 

where Ui are the displacement components: Ui := Xi — Xi. Here x is the 
current position of a particle initially at X. 



The equations of motion, in the absence of body forces, are given by 

dijkid'^Uk/dxidxj = pd'^Ui/dt^, (2.4) 

where p is the mass density of the crystal. 
We consider displacements of the form 

u = { Ae*'"(^"-*)}+ = e-"^"-"{ A+ cos u;(S+-x - t) - A^ sin u;(S+-x - t) }, 

(2.5) 
where S = S"*" + zS is the slowness bivector (complex vector) [U |3] and 
A = A+ + iA~ is the amplitude bivector (complex vector). When A is 
parallel to S, that is when A = aS, where a is some complex number, the 
inhomogeneous wave is said to be "longitudinal" [SI [6]. When A is isotropic, 
that is when A ■ A = 0, the wave is circularly polarized, the ellipse associated 
with the bivector A being a circle [3]. 

Circularly polarized longitudinal inhomogeneous waves are those waves 
for which A and S are both isotropic and parallel. For such waves, 

S ■ S = 0, (2.6) 

or equivalently, 

S+.S+ = S--S-, S+-S- = 0. (2.7) 

Thus, the planes of constant phase are orthogonal to the planes of constant 
amplitude; the waves propagate in the direction of S"*" whilst the amplitude 
decays in the direction of S~. The particle paths are circles in the plane 
of A"*" and A~, or equivalently, because A = aS, in the plane of S"*" and 
S~. The sense of description of the circle is from S"*" towards S~, retrograde, 
similar to the sense of Rayleigh waves propagating close to the free surface 
of a semi-infinite isotropic elastic material. 

Inserting 02. 5p into (12. 4 p gives the propagation condition, 

Qik{S)Ak = pAi, Qik{S) = dijkiSjSi, (2.8) 

where Qik{S) is the acoustical tensor corresponding to the slowness bivector 
S. A systematic procedure for obtaining all solutions A, S of (12. 8 p has been 
introduced by Hayes [1] and is called the "directional ellipse method" or 
"DE-method" . It consists in writing S as 

S = NC, (2.9) 



where A^ is a complex number and C is any bivector of the form C = mm+in, 
with m, fi, two unit orthogonal vectors and m > 1. We call A^ the complex 
scalar slowness and C the directional bivector. Inserting (12.91) into (12.81) 
yields the eigenvalue problem, 

Q^kiC)Ak = pN-''A, (2.10) 

for the complex symmetric tensor Qik{C)- All solutions of the propagation 
condition may then be obtained by prescribing C arbitrarily and solving 
(I2.10p for N'"^ and A, which gives the slowness bivector through (12. 9p . and 
the amplitude bivector A (up to a complex scalar factor). When, for some 
eigenvalue N~'^, the corresponding eigenbivector A is parallel to the direc- 
tional bivector C, the corresponding inhomogeneous plane wave is said to be 
longitudinal. Thus, for longitudinal inhomogeneous plane waves the propa- 
gation condition becomes 

QikiC)Ck = pN~^C„ (2.11) 

for some A^~^. Hence, these waves are only possible for directional bivectors 
C such that Q(C)C is parallel to C, or equivalently, 

CxQ(C)C = 0, (2.12) 

or 

Qik{C)Ck _ Q2k{C)Ck _ Q3k{C)Ck , s 

T; — T; — 7^ ■ {1.16} 

Oi O2 O3 

Here we seek particular classes of materials which are such that longitu- 
dinal inhomogeneous plane waves are possible for any choice of C satisfying 
C ■ C = (so that C is of the form C = m + in). That is, we wish to deter- 
mine under which conditions on the stiffnesses dijki are CPLIPWs possible 
for all choices of the plane of C, or equivalently of the normal a = m x fi. 



3 Crystals for which longitudinal circularly 
polarized waves are possible in all planes 

3.1 Necessary and sufficient conditions 

Here we seek under which conditions on the stiffnesses dijku we have: C x 
Q(C)C = for all C = m+m, or equivalently, for all C satisfying C ■ C = 0. 



For convenience we adopt the Voigt [7] contracted notation for the elastic 
stiffnesses, 

di2 = (iii22, (iss = c^ssaa) c^45 = ^^2313; dm = diiiii etc. (3-1) 

With these notations, 

Qii(C) = diiC^ + ^666*2 + d^^C-^ + 2diQCiC2 + 2di^CiC^ + 2^566*2(73, 

Q22(C) = d^^Ci + (i22C*2 + '^44C'3 + 2(^26^*1 C*2 + 2(i64C*lC*3 + 2(i24C*2C*3; 
Q33(C) = (issC;^ + ^446*2 + '^33^*3 + 2(i45CiC2 + 2d-i^CiC^ + 2^34(726*3, 

(3.2) 
and 

Ql2(C) = (iieC*]^ + d2QC2 + ^45^*3 

+ (C?12 + 46)C*iC*2 + (t/l4 + rf56)C*iC*3 + (^25 + C?46)C*2C*3, 
Q23(C) = d^gCj^ + (i24C*2 + C?34C*3 

+ (C?25 + rf46)C*iC*2 + (46 + 45)C*iC*3 + (^23 + (^44)6*26*3, 

Q3i(C) = di^Ci + (i46C*2 + ds^C^ 

+ {du + 46)C*iC*2 + (t/i3 + 45)C*iC*3 + (^36 + C?45)C*2C*3. (3.3) 

Consider first C = (0, l,i). It is isotropic. For this C, conditions (12.131) 
become QikCk = and iQ2kCk = QskCk, which read, exphcitly, 

(^36 + 2^45 - c/26) - i{d25 

4(^34 - ^24) + i{d22 + C?33 - 2^23 " 4^44) = 0. (3.4) 

Hence, 

0^26 = '^36 + 2(^45, (^35 = (^25 + 2^46, (^24 = '^34, 4(^44 = ^22 + d^s — 2^23- 

(3.5) 

We then consider in turn C = (^,0,1) and C = (l,i,0). These choices 

yield conditions of the type (13. 5p . which may be read off from (13.51) on cycling 

the indices 1— )-2— )-3— )-l,4— )-5— )-6— >-4. The complete set of conditions 

obtained in this way is 

dlG = 0^26 = '^36 + 2(i45, '^35 = ^^15 = ^^25 + 2(^46, ^24 = (^34 = C^14 + 2^56, 
4ci44 = Ci22 + ^3 - 2^23, 4^5 = ^33 + du - 2c?i3, 4^66 = du + d22 - 2rfi2. 

(3.6) 



he - 


- C?35) = 0, 


i23- 


- 4^44) = 0. 


dsA, 


4(^44 = (i; 



We refer to materials whose stiffnesses satisfy these conditions as "special". 
It follows that "special" materials have at most twelve independent elastic 
stiffnesses. For instance, they are 



dii di2 (ii3 (ii4 (ii5 diQ 

d22 d23 d24 C?25 • 

dss • • d-sQ 



(3.7) 



d34 


— d24, 


45 


= di5, 


d26 


= die, 



and the remaining elastic constants (denoted by "•" above) are determined 
from those 12 as 

dM = i{d22 + c^ss " 2^23)) '^45 = 2 ('^26 ~ c^se)) 

'^55 = lydii + (iss — 2(ii3), d^Q = 2 ('^34 — c^u) 

dm = i{dii + d22 — 2(ii2), d^^ = 2 ('^35 ~ '^25)- (3-8) 

We note that out of the nine conditions 03. 6p . six are "structurally invariant" 
mj for some rotations of the coordinate system. In particular, if the two 
conditions 

die - C?26 = dii + ^22 - 2rfi2 - 4^66 = 0, (3.9) 

are satisfied in the coordinate system linked to the crystallographic axes 
(0x10:2X3), then they are also satisfied by the stiffnesses d*j obtained from 
the dij after any rotation of the coordinate system about the X3 axis. These 
invariants are Type IB in Ting's classification [8]. Similarly, the conditions 

c^i5 - ^5 = dii + 43 - 2rfi3 - 445 = 0, (3.10) 

are invariants under rotation of the coordinate system about the X2 axis, and 

^4 - ^4 = d22 + 43 - 243 - 444 = 0, (3.11) 

are invariants after rotation of the coordinate system about the Xi axis. 

When the relations (13. 6p hold, it may be checked by direct calculation, 
using (13. 2p . (13. 3p . and taking into account the relation Cl + C2 + C^ = 0, 
that 

Q,fc(C)Cfc = piV-2(C)a, (3.12) 



for all isotropic bivectors C, with N^^(C) given by 

pN^\C) = lidnCl + ^22^2^ + rf33C|) + 2(rfi6CiC2 + rfssCiCa + ^24^2^3). 

(3.13) 
It follows that the relations (13. 6 p between the elastic stiffnesses are the nec- 
essary and sufficient conditions for CPLIPWs to propagate for all choices of 
isotropic directional bivectors, or, equivalently, for all choices of the polar- 
ization plane. 

The expression (I3.13P may also be written 

pN£\C) = C-^lC, (3.14) 

where 'I'l is given by 



*; 



die 2*^22 d24: 
d35 d24 2'^33_ 



(3.15) 



Associated with ^l is the "^L-ellipsoid" : x ■ (^^ + pl)x = 1, where p is 
chosen such that (^l +pl) is positive definite. If C is chosen to lie on either 
plane of central circular section of the ^L-ellipsoid, then C • ^lC = and 
thus N£'^(C) = 0: the corresponding waves do not propagate. 

In the next section, we derive the general structure of the acoustical tensor 
for "special" materials. 

3.2 Acoustical tensor 

From the previous section, we know that the acoustical tensor Qik{C), with 
C • C = 0, of "special" materials admits pA^^^(C) given by (I3.13P as an 
eigenvalue. The corresponding eigenvector C is isotropic, and so [1] this 
eigenvalue is a double eigenvalue. It follows that the other (simple) eigenvalue 
piV,"2(C), say, is given by pN^^C) = tr Q(C) - 2piV-2(C), that is, 

pN~^{C) = (t/55 + 46)^1^ + (C?44 + 46)^2^ + (C?44 + c/ss)^! 

+ 2d45ClC2 + 2d46ClC3 + 2rf56C2C3. (3.16) 

The expression (13.160 may also be written 

pN;\C) = C-^tC, (3.17) 
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where "^t is given by 



*i 



d55 + d^Q 
dm 



dm dm 

(^44 + dQQ ^56 

^56 C?44 + C^SS 



(3.18) 



Associated with ^t is the "^^^-elhpsoid" : x ■ {^t + pl)x = 1, where p is 
chosen such that ('i'T + pl) is positive definite. If C is chosen to he on either 
plane of central circular section of the ^T-ehipsoid, then C ■ ^^-C = and 
so A^~^(C) = 0: there is no corresponding transverse circularly polarized 
propagating wave. 

Now we compute the components of the matrix r(C) := Q(C)— pA^~^(C)l. 
We find, using the conditions (13. 6 p and Cf + Cl + Cl = 0, that 



(3.19) 





rn(C) = (^ + z/ - \)Cl + 27C1C2 + 2/3C1C3, 




r22(C) = {v + \- fi)Cl + 2aC2C3 + 27C1C2, 




r33(C) = (A + /i - u)Cl + 2/3C1C3 + 2aC2C3, 


and that 





where 

A: = 

a :- 



ri2(C) = -7CI + Z/C1C2 + aCCs + /3C2C3, 
123 (C) = -aC^i + AC2C3 + (3C1C2 + iCiC-i, 
131 (C) = -pel + /iCiC3 + 7C2C3 + «CiC2, 



2 ('^22 
(ii4 + 



d 



dss 

56) 



2rf. 



2 ('^11 + '^3 



44;, /^ :— 2l"ll "•" "33 
/3 := (^25 + '^46) 



2d, 



55 J 



Let M be the real symmetric matrix defined as 

A 



M: = 



-7 
-7 fi 
—(3 —a 



-/3 

—a 

z/ 



It can be checked that 



7 := 



2\dii 



2 

'^36 



(3.20) 



+ ^22 - 246) 

dm- (3.21) 



(3.22) 



r(C) = (A + /i + z/)C®C-MC®C-C® MC. 



(3.23) 



Alternatively, introducing M := |(A + yU + z/)l — M, the matrix r(C) may 
be written as 

r(C) = MC ® C + C ® MC, (3.24) 

where, explicitly, 

2(^11 + ^^44 — C^55 — f^ee) 1^36 + ^45 C?25 + '^46 

M= d3e + d45 ^{d22 + d55 - di4 - dm) du + d^ 

d25 + d46 di4 + dse 2 ("^33 + dsQ — du — d^5)_ 

(3.25) 
Hence, noting that C ■ C = and recalling the definition of r(C), we find 
that the acoustical tensor may be put in the form, 

Q(C) =piV,^2(C)l + MC®C + C®MC. (3.26) 

This decomposition of the acoustical tensor shows directly that the isotropic 
bivector C is an eigenvector of Q(C), whose eigenvalue pN^'^{C) given by 
( ]3.13p . can equivalently be written 

pAr-2(c) = pN;\C) + C ■ MC. (3.27) 

Also, the eigenvector corresponding to the eigenvalue pN~^{C) is C x MC. 
We call the corresponding wave the "transverse" inhomogeneous plane wave, 
because its amplitude bivector A is orthogonal to C: C ■ A = 0. In general it 
is elliptically polarized because (C x MC) ■ (C x MC) = -(C ■ MC)^ ^ 0. 
Of course C ■ [C x MC] = 0, which means |1] that the orthogonal projection, 
upon the plane of C, of the ellipse associated with the amplitude bivector 
C X MC, is a circle. 

Let C be a choice of C for which C ■ MC = C ■ C = 0, so that C is 
parallel to C x MC and there is a triple eigenvalue pN£^{C) = pN~'^{C). 
This special case occurs when the plane of C is one of the two planes of 
central circular section of the the "M-ellipsoid" : x ■ (M + pl)x = 1, where p 
is chosen such that (M + pi) is positive definite. 

Remark: A relationship between the scalar slownesses of certain 
waves 

The form 03.271) of the relation triggers the following remark. Let the matrix 
M have eigenvalues pi, P2, Ps and corresponding unit orthogonal eigenvectors 
ai, a2, as so that 

M = piai (g) ai + p2Si2 ® as + paag ® ag, a.^ ■ a^ = Sij. (3.28) 
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Then consider the following isotropic bivectors Ci, C2, C3, 

Ci:=a2 + za3, C2:=a3 + zai, C3 := ai + za2. (3.29) 

Clearly, 

CrMCi=p2-P3, C2-MC2=P3-Pi, C3-MC3=Pi-p2, (3.30) 

so that 

3 

J2Ci-MC, = 0. (3.31) 

1=1 
Then, from (I3.27P and (I3.3ip . we have the relation 

3 3 

E^^~'(CO = E^*~'(CO' (3-32) 

1=1 i=l 

for the Cj given by (K^ . 

Example 

To illustrate the results of this section, we work out a simple example. Let 
C = i + ij. Then the corresponding acoustical tensor Q(i + ij) is 

c^ii - ^66 + 2idiQ i{dQe + ^12) di5 - d^Q + i{d^Q + du) 

i(46 + c^i2) dee - d22 + 2id2e -{d56 + du) +i{di5 - d^e) 

di5 - diQ + i{d5Q + du) -id5Q + du)+i{di5-d4Q) d^^ - d^^ + lid^rj 

(3.33) 
Computing Q(C)C and using the conditions (13. 6p . we find that C is indeed 
an eigenvector of the acoustical tensor, with eigenvalue given by (I3.13p . 

Q(C)C = pN-'C, pNf = |(rfn - ^^22) + 2td,,. (3.34) 

Further, computing C x MC and then Q(C)(C x MC) and using the con- 
ditions (13. 6p . we find that C x MC is also an eigenvector of the acoustical 
tensor, with eigenvalue now given by (13.160 . 

Q(C)(C X MC) = piV-2(C X MC), pN;^ = rfgs - ^4 + 2id45. (3.35) 

In this simple example, we prescribed the normal to the plane of the 
isotropic slowness bivector (or equivalently, to the plane of the amplitude 
bivector) to be k. Then, we chose C to be i + zj, leading to a complex 
eigenvalue. In the next section, we show that it is always possible to choose 
C such that the corresponding eigenvalue is a real positive number. 
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Remark: General form of the acoustical tensor for the "special" 
materials 



Using the relations fl3.6p in the expressions Qij{C) given by (13. 2 p and (13. Sp . it 
may be seen that without any restrictions on C, the acoustical tensor Q(C) 
for the "special" materials may be written as 



Q(C) = pN;\C)l + MC ® C + C ® MC + (C ■ C)A. 



(3.36) 



Here the expression for pN^ ^(C) is given by (13.161) and A is defined by 



— Ci44 045 U4g 

dm d^Q —d^e 



Comparing (I3.18P and (13.371) . we note that ^t = A 

bydXIZD, 

pN;\C) = C ■ AC - (tr A)C ■ C. 



(3.37) 

(tr A)l, so that 
(3.38) 



It follows that 

Q(C) = [C- AC-(tr A)C-C]1 + MC®C + C®MC + (C-C)A, (3.39) 

which is an expression for the acoustical tensor Q(C) for "special" materials, 
given in terms of only two matrices, A and M. 

We note for "special" materials that dijki may be written 

dijkl = ^jl^ik + \kSjl + \l^jk + ^jk^il — ^ij^kl — ^kl^ij 

- (tr A){6ikSji + 6jkSii - Sij6ki) + MijSki + MkAj- (3.40) 

Remark: Homogeneous plane waves in "special" materials 

For homogeneous plane waves propagating in the direction n in the "special" 
materials, let u = Aexp2a;(S'n ■ x — t). The corresponding propagation 
condition is then 

Q(n)A = pS^A, (3.41) 

where Q(n) is given by (I3.39P with ni, n2, n^ replacing Ci, C2, C3, respec- 
tively. Thus 



Q(n) = [n ■ An - (tr A)]l + Mn (g) n + n (g) Mn + A. 



(3.42) 



12 



We note a special solution. 

Let the Hamiltonian decomposition of A be given by [3] 

A = /il + K(h+®h" + h~(8)h+), (3.43) 

where /x, k are constants and h^ are constant unit vectors. Choose n = n* 
such that 

(n* X Mn*) X (h+ x h") = 0. (3.44) 

(It is always possible to do this.) A solution for a homogeneous plane wave 
propagating along n* is 

u = (h+xh")expza;(A^"^n* ■ x-t), where piV~^ = /x+n* ■ An*- (tr A). 

(3.45) 

4 Description of the waves 

Before proceeding, we note that the two equations (I3.13P and fl3.16p . giving 
the complex scalar slowness of the longitudinal circularly polarized wave and 
of the transverse elliptically polarized wave, respectively, have the same form, 

C-^C = pN-\ C-C = 0, (4.1) 

where ^ is a real symmetric tensor given by ^ = ^l (see (13.151) ) for the 
longitudinal wave and by ^ = ^t (see (I3.18P ) for the transverse wave. 
Accordingly, we determine the details of the slowness of the corresponding 
wave solutions. The results may be adapted either to the longitudinal wave 
or to the transverse wave by replacing ^ with either ^l or ^t- 

4.1 Construction of the slowness bivector 

Let \l/i, \l/2, "^3 be the eigenvalues of ^ and ei, 62, 63 the corresponding 
orthogonal unit eigenvectors. We assume that the eigenvalues are ordered as 
\l/i > \l/2 > ^3- We note that the pair (14.11) is equivalent to the pair 

C • (* + pl)C = pN-\ C ■ C = 0, (4.2) 

where p is an arbitrary constant. By choosing p suitably large and positive, 
so that \E'3 +p > 0, we may define the positive definite matrix (^ +pi) and 
associate the ellipsoid 

x-(*+pl)x = l, (4.3) 
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with ^. We call this the "^-ellipsoid" . The planes of central circular section 
of the ^-ellipsoid have unit normals h^ given by [3] 

i , i 1 

(^1 - ^3) 2 h± = (^1 - ^2) 2 ei ± (^2 - ^3) 2 eg. (4.4) 

We call these normals h^, the "optic axes" of the ^-ellipsoid and note that 
they are independent of p. 

To describe the slownesses of the waves corresponding to (14. ip . we recall 
that 

S = A^C, C = m + in, mm = nn = l, mn = 0. (4-5) 

Here C is such that C ■ C = and thus S ■ S = 0. Also (m, n) is any pair of 
orthogonal unit vectors in the plane of S, or equivalently in the plane with 
unit normal a := m x n. In (14. 5p . A^ is a scalar to be determined from the 
equations (14. ip . 

Because we are at liberty to choose for (m, fi) any orthogonal unit pair 
in the plane with normal a, we choose (m, n) along the principal axes of the 
elliptical section of the ^-ellipsoid by the central plane a ■ x = 0. Specifically, 
we take rii as the unit vector along the minor axis of this ellipse and n along 
the major axis. In that event, 

m ■ *n = 0, (4.6) 

and using (14. 5 p and (14. 6p . the pair (14. ip give 

pN'"^ = m ■ *m - n ■ *n > 0, (4.7) 

in general, so that A^~^ is purely real: A^~^ = v, say. Also, 

S+ = v-^ui, S~ = v~^n, (4.8) 

so that the planes of constant phase (amplitude) are: m ■ x = constant 
(fi ■ x = constant). 

Of course, if m • ^m = n ■ ^n, so that the radii to the ^-ellipsoid along 
the orthogonal unit vectors are equal, then the plane a ■ x = is a plane of 
central circular section of the ^-ellipsoid, and from (14.71) there is no propa- 
gating solution: pN'"^ = 0. 

Thus, in general, the slowness bivectors corresponding to (14. ip are ob- 
tained by first determining the central elliptical section of the ^-ellipsoid by 
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the plane a ■ x = 0. Then m and fi are chosen along the principal axes of 
the ellipse, and S+ and S~ are given by fl4.7p f l4.8p . 

To complete the picture we recall the results |3] for the determination of 
the principal axes of the central elliptical section of the ^-ellipsoid by the 
plane a ■ x = 0. 

In the determination there are three cases to be considered: 

Case(i) Normal a not coplanar with the optic axes. 

Case(ii) Normal a coplanar with the optic axes but not parallel to either 
optic axis. 

Case(iii) Normal a parallel to an optic axis. 

We can dispose of Case (iii) immediately because we have just seen that 
there is no propagating circularly polarized solution when a ■ x = is a plane 
of central circular section of the ^-ellipsoid. 

In Case (i), a ■ h+ x h~ 7^ 0, it has been shown [3] that m, n are in the 
direction of r^ given by 

r± = [h+ - (cos(/)+)a]/(sin(/)+) ± [h^ - (cos 0") a] /(sin (/)"), (4.9) 

where 0^ is the angle between a and the optic axis h^. Essentially, m and 
h are along the internal and external bisectors of the angle between the 
orthogonal projections of h+ and h" onto the plane a ■ x = (this is the 
Fresnel construction of Optics.) 

In Case (ii), a ■ h+ x h~ = 0, a x h^ 7^ 0, it has been shown [3] that 

m = [h+ - (cos0+)a]/(sin0+), fi = (h+ x a)/(sin(/)+). (4.10) 

Essentially, m is along the orthogonal projection of h+ (or h~) onto the plane 
a ■ x = 0, and n is orthogonal to that plane. 

This means that as the direction of the unit normal a is varied, the 
corresponding m and n are given by (14. 9 p and (I4.10p . Also, as shown by 
Boulanger and Hayes [3], in Cases (i) and (ii), m ■ ^m — fi ■ ^n = (\E'i — 
\l/3) sin0+ sin0~, so that from (14. 7p . 

pN""^ = pv"^ = (^1 -^3)sin0+sin0", (4.11) 

where 0^ are the angles that the normal a to the plane of S makes with the 
optic axes. 

We summarize the situation. 

To determine all S"*" and S~ corresponding to (14. ip . the "optics axes" h''^ 
are first determined. Then the plane of C = rh + m is chosen; it has unit 
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normal a. If a is not coplanar with the two optic axes, then the corresponding 
m and fi are in the directions of r^ given by f l4.9p and A^ is given by fl4.7p . If 
a is coplanar with the optic axes h^ but not along either of them, then the 
corresponding m and n are given by f l4.10p and A^ is given by (14. lip . If a is 
along h+ or h~, there is no propagating wave. 



4.2 Example: Crystal with a plane of symmetry 

To illustrate the method described above, we take a "special" material with 
a symmetry plane at X3 = 0, say. Then the tensor ^l defined in ( I3.15P for 
the longitudinal wave is given by 



*; 



kdu 


die 





dw 


^,d,2 











y 



33 



(4.12) 



Its eigenvalues are 



^1, 



du + d22 ± \ {dn - 0^22)^ + 16rf 



v]/. 



odsa- 



(4.13) 



Here we assume that the stiffnesses of the material are such that these eigen- 
values are ordered \E'i > ^2 > ^3- The corresponding unit eigenvectors 
are 



"|C?22 - ^1" 




"0" 




1 

'6 


-die 


-die 


, 62 = 





, 63 = 


\dii - *3 







1 








^' = -5 



where 5 is the positive quantity given by 



(4.14) 



5' 



\{dii - d22f + 2dl^ + \{dii - d22)^{dii-d22? + lQdl^. (4.15) 



With this choice, the optic axes defined by (14. 4 p lie in the symmetry plane 

X3 = 0. 

For simplicity, we now focus on waves polarized in the symmetry plane, 
that is, we choose the normal a to the plane of S to be along 62. Then 
a is obviously not coplanar with the optic axes (Case (i) of the previous 
subsection). Specifically, the angles between a and the optic axes are (j)'^ = 



16 



7r/2. It follows from fl4.9p that m and fi are in the directions of h+±h 



I.e. 



m = ei, n = 63. (4.16) 

We conclude that the following CPLIPW may propagate in a monoclinic 
crystal with symmetry plane at X3 = and with stiffnesses satisfying fl3.6p . 

u = e"''''^'''{ei cos /c(ei-x - vt) - 63 sin /c(ei-x - vt)}, (4-17) 

where the orthogonal unit vectors ei, 63 are defined in fl4.14p . k is an arbitrary 
real wave number, and the real speed v is given by 



pv^ = ^1 - ^3 = l^idu - d22Y + 16rf?6. (4.18) 

5 Crystals with symmetries 

In this section we investigate how the conditions (13. 6p for an anisotropic 
crystal to admit a CPLIPW for all choices of the polarization plane are 
affected when the crystal presents certain symmetries. 

5.1 Monoclinic crystals 

Here we consider crystals with a plane of symmetry, at X3 = say. For that 
class of materials, 

(ii4 = (ii5 = (i24 = (^25 = '^34 = '^35 = (^46 = ^^56 = 0. (5-1) 

It follows that four out of the nine equations (13. 6p reduce to trivial identities, 
automatically satisfied. The nine conditions (13. 6 p reduce to a set of five 
equations, 

dl% = (^26 = C?36 + 2(i45, 

Mm = d22 + d-is - 2^23, 4c?55 = ^33 + dn - 2di3, 4c?66 = du + ^22 - 2di2. 

(5.2) 

5.2 Orthorhombic crystals 

Orthorhombic crystals possess three symmetry planes, at Xi = 0, X2 = 0, 
X3 = 0. In addition to (15. ip . the relations 

die = d26 = dsQ = di^ = 0, (5-3) 
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also hold. The set of nine conditions 03. 6 p now reduces to a set of three 
equations, 

4^44 = d22+d33-2d23, 4^5 = (igg +(iii - 2(ii3, 4^6 = ^11 + 1^22 " 2^12- (5.4) 

5.3 Trigonal, tetragonal, and cubic crystals 

For trigonal crystals, ^24 = —du = — c?56 7^ 0, and one of the conditions ( 13. 6p . 
namely: ^24 = (ii4 + 2^55, cannot be satisfied. 

For tetragonal crystals, dn = ^22, c?i3 = c?23; (^44 = c^ss? and the condition: 
4(i66 = c^ii + '^22 — 2(ii2 reduces to: d^Q = {dn — di2)/2, which would mean 
that the crystal is in fact hexagonal (transversally isotropic). 

For cubic crystals, dn = ^22 = c?33, di2 = (^23 = '^13, (^44 = c^ss = d^Q, and 
the condition: 4^66 = dn + 6/22 — 2(ii2 reduces to: (ige = (c^ii ~ '^i2)/2, which 
would mean that the material is in fact isotropic. 

We conclude that there are no trigonal, no tetragonal, and no cubic crys- 
tals in which CPLIPWs may propagate for all orientations of the slowness 
plane. 

5.4 Hexagonal crystals 

For hexagonal crystals, the following relations hold for the stiffnesses, 

du = d22i di3 = d23, (^44 = (^55, (^66 = (^^11 ~ di2)/2, (5.5) 

in addition to (15.1 p and f l5.3p . The nine equations ( 13. 6 p reduce to a single 
equation, 

du = iidn + d33-2di3). (5.6) 

5.5 Isotropic materials 

For isotropic materials, the following relations hold for the stiffnesses, 

du = d22 = d33, di2 = d23 = di3, (^44 = ^55 = d^Q = [dii — di2)/2, (5.7) 

in addition to ( 15. ip and ( 15. 3p . Then, the nine equations (13. 6p are all identi- 
cally satisfied. However, the propagation condition (I3.13p . giving the complex 
scalar slowness Nl now simplifies to 

pN^'iC) = idn(Ci2 + C, + C|) = 0, (5.8) 



for isotropic slownesses. It follows that CPLIPWs may not propagate in an 
isotropic material, for any choice of isotropic slowness. 

On the other hand, this analysis shows that "longitudinal" static expo- 
nential solutions with an isotropic slowness bivector always exist for linear 
isotropic elastic materials. This result may be checked directly in the fol- 
lowing manner. Recall the classical equations of equilibrium of an isotropic 
material, 

(A + fi)uj^ij + fiUijj = 0, (5.9) 

where A and /i are the Lame constants. Then it is easy to check that the 
field 

«. = ^.e-^.-., S,S, = 0, (5.10) 

is indeed an exact solution. 

We sum up the situation. There are no isotropic, cubic, trigonal, or 
tetragonal elastic crystals such that CPLIPWs may propagate in every plane. 
However, the propagation of CPLIPWs is theoretically possible for all planes 
of some triclinic, monoclinic, or orthorhombic elastic crystals provided some 
relations among the elastic stiffnesses are satisfied. We note that there is 
unfortunately insufficient data for such crystals available at present to enable 
us present an explicit example, but we recall that the values of the elastic 
stiffnesses change with pressure, temperature, prestress, etc. and that they 
may consequently be adjusted to produce an adequate crystal. 
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